Fano resonance in a two-level quantum dot side-coupled to leads 
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We theoretically study Fano resonance in a two-level quantum dot side-coupled to two leads, 
which are connected by a direct channel. The resonance lineshape is found to be deformed, from the 
conventional Fano form, by interlevel Coulomb interaction and interlevel interference. We derive 
the connection between the lineshape deformation and the interaction-induced nonmonotonicity of 
level occupation, which may be useful for experimental study. The dependence of the lineshape on 
the transmission of the direct channel and on the dot-lead coupling matrix elements is discussed. 
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Fano resonance,^ which is the interference between a 
resonant state and a continuum, appears ubiquitously in 
various systems. It has a lineshape of the form 

Here, e = (cj— ^q) /F is the detuning parameter measuring 
energy uj from the resonance center and normalized by 
the resonance half-width F, and q is the Fano parameter 
characterizing lineshape asymmetry. For g — > oo Eq. ([T]) 
becomes the Breit-Wigner form, while for q = it shows 
an antiresonance. In general, g is a complex quantity.— 

Recently, Fano resonance has been investigated in 
mesoscopic electron systems such as waveguides,^ quan- 
tum dots^i^ Aharonov-Bohm rings coupled to a quan- 
tum dot,^'^'^ and carbon nanotubesj^^ The studies im- 
ply that Fano resonance provides a useful tool study- 
ing dephasingi^ On the other hand, some aspects of 
the interplay between Fano resonance and Coulomb in- 
teraction have been studied. They include Fano reso- 
nance modified by charge sensin g^^d^ and Fano-Kondo 
antiresonance in a spin-degenerate single-level quantum 

The Fano lineshape ^ is applicable for a system with 
a single resonant level. It is valid as well for multilevel 
systems as long as each single-particle level is well sep- 
arated from the adjacent levels in energy. Most stud- 
ies on Fano resonance have been carried out mainly in 
this single-level regime. However, one may often find 
the multilevel regime where single-particle level spacing 
is comparable to level broadening. In this regime, the 
single-level Fano form is not applicable any more. More- 
over, this regime possesses interesting effects, absent in 
the single-level regime, such as nonmonotonic level oc- 



cupation 
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due to Coulomb repulsion. It has been 



reportedi^ that the nonmonotonic behavior of level oc- 
cupation influences Breit-Wigner lineshape. Therefore, it 
may be interesting to see the modification of resonance 
lineshape, from Eq. ([1]), in a more general multilevel Fano 
regime and to analyze the infiuence of Coulomb inter- 
action on resonance lineshape, which is the aim of the 
present paper. 

In this work, we theoretically study Fano resonance in 
a two-level electron quantum dot (QD) side-coupled to 



two leads, which are connected by a direct channel (see 
Fig. [T|) . Interlevel Coulomb repulsion in the dot is taken 
into account and the spin of electrons is neglected for sim- 
plicity. We use Keldysh formalism and a self-consistent 
Hartree-Fock (SCHF) approach, to obtain and to ana- 
lyze the Fano resonance lineshape of the two-level sys- 
tem. The two-level lineshape is found to be deformed 
from the single-level form by the Coulomb repulsion 
and interlevel interference. We derive the connection 
[Eqs. to between the Hneshape deformation 

and the nonmonotonicity of level occupation, when the 
QD level spacing (after renormalized by the repulsion) is 
larger than level broadening so that the nonmonotonic- 
ity is not too strong. The connection may be useful for 
experimental study. We also discuss the dependence of 
the nonmonotonicity on the direct-channel coupling, an 
extension to a spinful single-level case, and the tempera- 
ture range where the SCHF result is valid, below which 
correlation-induced resonance a^°i^^'^^i^^i^^ may emerge. 

We start with the Hamiltonian Hu of the spinless two- 
level QD, 



Hd = ^ (f^-eyg)4d^ + [/44d2di, (2) 

7=1,2 

where dt^ creates an electron at QD level 7 = 1,2 with 
energy Vg is the gate voltage apphed to the QD, 
and U is the interlevel Coulomb repulsion. Without 




FIG. 1: A quantum dot with spinless two levels 7 = 1,2, side- 
coupled, with coupling matrix element ti:,(_B)-y, to two leads L 
and R, which are connected by a direct channel with coupling 
tLR. There is Coulomb repulsion U between the two levels. 
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loss of generality, we can set < ^2- The QD cou- 
ples via tunneling to noninteracting leads a = L^R, 
which are connected by a direct channel (see Fig. [1]), 
so that the total Hamiltonian of the system is H — 
Hd + Hl + Htd + Htl, where the two leads, the dot-lead 
tunneling, and the lead-lead tunneling are described by 

and Htl = J2kk'i^LRcli^Ck'R + h..c.), respectively. Here 
^ka creates an electron with momentum k and energy £^ka 
at lead a. We ignore for simplicity the momentum de- 
pendence of tunneling matrix elements ta-y between level 
7 and lead a, and that of ti/j of the direct channel. 

It is worthwhile to note that electron transport through 
the two levels depends on the dot-lead coupling taj. To 
see this, we consider simple cases^ with time reversal 
symmetry where tfj^'s are chosen to be real, i^,! = tiu = 
ti > 0, = stR2 = ^2 > 0, and the phase parameter 
s = ±1. After redefining a pair of two orthogonal leads, 
saying a = ±1, by Ck& = {ckL + ackB)l^/2, one finds 
that for s = +1, the two QD levels couple to the same 
lead a = -fl, while for s = — 1, they couple to different 
leads each other. This s-dependent nature of coupling 
to leads a characterizesi^iiSiS^ electron transport such as 
interference. Below, we will use the above choice of tc^'s 
and see the dependence of the two-level Fano lineshape 
on s. Note that Ilr is also chosen to be real as well. 

We obtain electric current of the QD system using 
Keldysh formalism. The current, Jl — —e{fiL) — 
{ie/h){[nL, H]), in the lead L can be expressed as 

k'y 



kk' 



(3) 

Here, = c^i^CkL is the electron density operator, 
and G^(u;)'s are lesser Green's functions which corre- 
spond in time domain to ^.j^{t,t') = i{c\j^{t')d^{t)) 

and G^,^ kLi^'^') ~ *(Cfc_L(^')cfc'i?(i))- finds the cur- 
rent Jr in lead R in the same way. After some algebra3^ 
using the relations connecting G^'s and the retarded 
Green function G;^^, of the QD [see Eqs. (gl) and dH)], the 
wide-band approximation for lead Green's function gka, 
J2k 9ka i'^) — P being the density of states of leads, 

and the steady-state current conservation — —Jr^ we 
arrive at a useful form of the current, 



du;{fL{cj)-fR{cj)){TB+TQD{u;)), (4) 



where fa is the Fermi distribution of lead a. The back- 
ground transmission, Tb — An'^ p'^tj^j^/x'^ , comes only 
from the direct channel, where x = l+n^p^tj^j^ is the fac- 
tor counting multiple reflections via the direct channel. 
The term Tqd results from the paths passing through 
the QD, and it depends on the phase parameter s, 



'QD 



-am[z/fiGli +2/T2G^2 +z2yfif2GliG^*] , 

(6) 

where y = (1 — iirptLR)'^ / x"^ is a complex factor coming 
from the effect of the direct channel and 



f -^2- 



(7) 



is the QD level broadening i^^i^^ Notice that becomes 
narrower, as Ilr increases, from the level broadening 
F^ = 27rpi^ of the QD without the direct channel. The 
first two terms of Eqs. (|5l6p describe the direct contri- 
bution through the QD level 7 as well as the interference 
between the paths through the level 7 and the direct 
channel, while the third shows the interlevel interference 
between the paths through the two levels. 

The derivation oi Tqe, depends on s due to the cou- 
pling nature to the leads a. For s — ~\, GI2 does 
not appear in Tqd, a-nd it is necessary to use Keldysh 
equation G^ = G''I]^G° for the dot lesser function. We 
use the noninteracting form of the lesser self-energy E< 
of the dot coming from the lead-dot coupling, therefore 
'Tq~^~^ in Eq. ([6]) is an approximate form valid within 
the SCHF approach used below. On the other hand, for 
s = the Keldysh equation is not necessary in the 
derivation, thus ^5^^ in Eq. ^ is exact. Note that 
Eqs. (|4l5|6p are reduced into the forms found in the pre- 
vious works on the single-level Fano resonance^'^'l^ (when 
S.2 — £.1 3> F^=i.2) and on a multilevel QD without the 
direct channel^^ {t^R — > 0). 

We obtain the retarded QD Green's function G^y and 
the level occupation {n^-y') in equilibrium by using the 
equation of motion method and the SGHF approach. 



G^2(^) 

('^77') 



^{uj--^^ + eVg-U{n^) 
I dujf{uj)3mG''^,^{uj), 



(8) 
(9) 
(10) 



(5) 



where D = (cj - + eVg - [/(na) - (a; - 6 + eVg - 

t/(ni) — S2) — (f7(ni2) — 115^2)^' 7 means the level different 
from 7, {n~f) = {n^-,) is the occupation of level 7, and 
f{Lo) = l/{eJ^'^ + 1). The self-energies are found to be 
S!^ = -{s'^^^'KptLR + and YT^^ = SsA-^pthR + 
i)\/f if2; from UmEIJ,, one can get Eq. We remark 
that Y7i2 and {ni2), therefore G5^2i vanish for s = — 1, 
because of the s-dependent coupling nature to the leads 
a. The SCHF approach is good when F^ is not too large 
compared with level spacing ^2 — Ci We later discuss 
the temperature range where the SCHF result may be 
valid. Since we have interest in resonance lineshape, we 
will focus on the the linear response regime below. 

We first discuss two-level resonance lineshape in the 
noninteracting case of [/ = 0. The two-level lineshape 
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FIG. 2: (color online) Two-level Fano lineshape T, level oc- 
cupation (n-y), and occupation cross-correlation {Sni5n2) as a 
function of Vg for (a) s = — 1 and (b) s = -f 1 in the noninter- 
acting case of (7 = 0. We choose Fi/F = 0.63, F2/F = 0.37, 
(6-6)/r = 1.6, and TvptLR = 0.3, where F = F1+F2. Zero- 
temperature results of {n^^i) and T are used for simplicity. 



can be obtained as T = Tb + 7qd(ijJ = 0), 

- (C^ - l)5s,i + sq{ei + 362 - K5s 



T = Tp. 



\ei£2 



- (C2 - l)<5,,i)2 (£1 + S£2 - 2C<5.,l)' ' 

where e-^ — {eVg — ^7)7^7 + s'' ^irptLB. is the detuning 
parameter of level 7 = 1,2, ( = -(fif2)"i/2fj^eS^2 = 
irptLRSs,!, the terms with Kronecker delta Sg^i come from 
S5^2, q = ^'(l — Tb)/Tb, s'^~^q is the Fano parameter of 
the level 7. Note that the lineshape T is reduced into the 
single- level form ^ when the level spacing £,2 — ^1 *C f^. 

Figure [2] shows typical two-level Fano lineshapes as a 
function of gate voltage Vg for the cases where U = and 
the level broadening F^ is comparable to the level spac- 
ing £2 — The entire lineshape may be understood as 
s-dependent mixture of interferences between the paths 
through one QD level and the direct channel (character- 
ized by the Fano parameter of the level) and between the 
paths through the two QD levels. For s = — 1, the upper 
resonance at ^2 has a negative value of Fano parameter, 
while the lower one at ^1 has a positive value. Therefore, 
the two resonances are out of phase (with difference by 
tt), as shown in Fig.[2fa). On the other hand, for s = +1, 
the two resonances are in phase [see Fig. [2Kb)]. 

In Fig. [51 we plot the cross-correlation of level oc- 
cupation,^''' {SniSn2) = (nin2) — (ni)(n2) ~ — (n-12)^, 
which may give more understanding of the s-dependent 
features. For s = +1, {Snidn2) can have a finite value 
and give rise to nonmonotonic behavior of (n-y) [see (ni) 
around eVg = ^2 in Fig. IHb)], while it vanishes for 
s — ~1. Note that it becomes suppressed as tLR in- 
creases since the overlap between the two levels or the 
level broadening F^ is reduced [see Eq. 0]. 

Hereafter we turn on the Coulomb repulsion U and 
discuss how it modifies the resonance lineshape. We find 
that within the SCHF treatment, the lineshape has the 
same form as the noninteracting case oi U = 0, Eq. pT|) . 



[£i£2 - s - ie - l)Sss + sqjii + SS2 ~ 2C<5,,i)]^ 

(£i£2 - s - iC - l)Ss.ir + (ffi + - 2C<5.,i)2 

(12) 
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FIG. 3: (color online) Upper panels: The same plots as in 
Fig. m in the interacting case of fZ/F — 25. Lower panels: 
The second resonances in (a) and (b) are analyzed using the 
noninteracting Fano form To and the nonmonotonicity A for 
(c) s = — 1 and (d) s = -1-1, respectively (see text). 



but with mean-field shifts 

. _ U 
Sj — Sj — ^ 



U 



{ni2 



(13) 



Note that Sj and (n^^/) depend on Vg. The shift of the 
detuning parameter can be understood as the level spac- 
ing renormalization due to the Hartree repulsion. On the 
other hand, the shift in ( comes from the Fock exchange, 
which is absent in the case of s = —1. 

In Figs. [Sja) and (b), we plot the lineshape T when 
nonmonotonic behavior— liiii^ii^ of (rij) occurs (see, e.g., 
{ni) around the second resonance). For s = -t-1, the 
nonmonotonic behavior comes from the Hartree repulsion 
as well as from the Fock exchange, while for s = —1, it 
is caused only by the former. Therefore, the s = -1-1 case 
shows the nonmonotonicity in a wider range of Vg where 
{SniSn2) is enhanced by the Fock exchange. 

The nonmonotonic dependence of {n-yj') on Vg modifies 
the lineshape from the noninteracting cases. Such mod- 
ification can be analyzed when the level spacing renor- 
malized by the Hartree contribution is much larger than 

U > f^, so 



level broadening F-,,, i.e., when £2 — £, 



that the nonmonotonicity is not too strong. In this 
case, the lineshape (|12p can be simplified, for gate volt- 
age, for example, around the second resonance (eVg ~ 
£,2 + U + sTrptLfiT2), into the single-level Fano form Tq, 



Tie2,q)^%{e2,q)^TB^-^"^^' 



1 



(14) 



but with not simple dependence of (ni) on Vg [see 



Eq. (uni)] 

measure 

MVg. 



To analyze further, we define nonmonotonicity 
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FIG. 4: (color online) Nonmonotonicity A of level occupation 
as a function of Vg for (a) s = — 1 and (b) s = +1 in the 
interacting case of f//r = 25. Different values of Ilr are 
chosen. The other parameters are the same as in Fig. |31 



Here, (711)0 (Vg) is the level occupation in the absence 
of the Coulomb repulsion, which can be obtained from 
Eqs. (|8l9ll0p . We can take {ni)o{Vg) ~ 1 in this case of 
large renormalized level spacing. For |A| <C 1, one has 
an approximated form of the lineshape, 



T 



^o(£2,o,9) 



d£2.0 



(16) 



The leading term Tq obeys the single-level Fano form ([T]) 
with the detuning parameter £2,0 = (eV^ — £,2 — 
U{ni)o)/f2 + sTvptLR ^ {eVg - - ^^)/f 2 + snptLR, 
which approximately linearly depends on Vg, and the sec- 
ond deformation term is proportional to A(V^), therefore 
it provides the connection between the lineshape and the 
nonmonotonicity. Around the first resonance, one can 
find the same forms of the connection and the nonmono- 
tonicity measure but with level index exchange 1 <-> 2 
and (712)0 — 0- The connection is applicable as well in 
the Breit-Wigner case of Tbw,o = %{tLR ^ 0)- 

In Figs. EJc) and (d), we plot the deviation of T from 
To as well as A. The deviation depends on the phase pa- 
rameter s, as the nonmonotonicity does. The lineshape 
deformation occurs in a wider range of Vg in the case of 
s — +1, since the Fock exchange enhances {dniSn2) and 



therefore additional nonmonotonicity only in the s = -1-1 
case, as discussed before. The connection between the 
lineshape deformation and the nonmonotonicity of level 
occupation found in Eq. p6p may suggest an experimen- 
tal study of the nonmonotonic behavior. 

We finally discuss a few remarks briefly. First, the non- 
monotonic behavior becomes weakened for larger direct 
channel coupling t^/j, because the level broadening F^ 
becomes narrower [see Fig. S] and Eq. ([7])]. Second, we 
extend the above findings to a spinful single-level QD. 
In this case, the Hartree repulsion induces the nonmono- 
tonic behavior of level occupation, but there is no inter- 
level interference and no Fock contribution. We find that 
the connection (|16p between the lineshape deformation 
and the nonmonotonicity is still hold (at temperatures 
larger than the Kondo temperature). 

Third, the temperature range where the SCHF ap- 
proach is valid depends on tLR. It has been foun d^^'^'^i^^ 
that when Ilr = 0, the two-level QD can be mapped into 
a Kondo system and it can show correlation-induced res- 
onance^S below the Kondo temperature. We find that 
the two-level QD with the direct channel can be also 
mapped^I, into a Kondo system, depending on s, and that 
the corresponding Kondo temperature decreases with in- 
creasing tLR. Therefore, our approach may be valid 
above the Kondo temperature, the upper bound of which 
can be estimated^^ from the case of t^u = 0. 

In summary, we have studied Fano resonance lineshape 
and the nonmonotonicity of level occupation in a two- 
level QD side-coupled two leads. The two-level lineshape 
is derived for both the noninteracting and interacting 
cases [Eqs. (|11I12I) ]. We especially obtain the connection, 
Eq. (|16p . between the nonmonotonicity and the Coulomb 
modification of the lineshape. We also find that stronger 
direct-channel coupling weakens the nonmonotonicity. 
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1 U. Fano, Phys. Rev. 124, 1866 (1961). 

^ A. A. Clerk, X. Waintal, and P. W. Brouwer, Phys. Rev. 

Lett. 86, 4636 (2001). 
^ E. Tekman and P. F. Bagwell, Phys. Rev. B 48, 2553 

(1993); J. U. Nockel and A. D. Stone, Phys. Rev. B 50, 

17415 (1994). 

J. Gores et al, Phys. Rev. B 62, 2188 (2000). 
^ I. G. Zacharia et ai, Phys. Rev. B 64, 155311 (2001). 
^ K. Kobayashi et al., Phys. Rev. Lett. 88, 256806 (2002). 
^ K. Kobayashi et al., Phys. Rev. B 70, 035319 (2004). 
® A. Fuhrer et ai, Phys. Rev. B 73, 205326 (2006). 
^ J. Kim et al, Phys. Rev. Lett. 90, 166403 (2003). 
-° W. Yi et ai, Phys. Rev. Lett. 91, 076801 (2003). 
^ A. C. Johnson et al, Phys. Rev. Lett. 93, 106803 (2004). 
^ P. Stefahski, A. Tagliacozzo, and B. R. Bulka, Solid State 

Commun. 135, 314 (2005). 
^ B. R. Bulka and P. Stefahski, Phys. Rev. Lett. 86, 5128 

(2001). 



W. Hofstetter, J. Konig, and H. Schoeller, Phys. Rev. Lett. 
87, 156803 (2001). 

M. Sato et al, Phys. Rev. Lett. 95, 066801 (2005). 

R. Berkovits, F. von Oppen, and Y. Gefen, Phys. Rev. 

Lett. 94, 076802 (2005). 

J. Konig and Y. Gefen, Phys. Rev. B 71, 201308 (2005). 
^® M. Sindel et al, Phys. Rev. B 72, 125316 (2005). 

M. Goldstein and R. Berkovits, New J. Phys. 9, 118 (2007). 
^° V. Meden and F. Marquardt, Phys. Rev. Lett. 96, 146801 

(2006) . 

C. Karrasch, T. Enss, and V. Meden, Phys. Rev. B 73, 
235337 (2006). 

^2 H.-W. Lee and S. Kim, Phys. Rev. Lett. 98, 186805 (2007). 
V. Kashcheyevs et al, Phys. Rev. B 75, 115313 (2007). 
P. G. Silvestrov and Y. Imry, Phys. Rev. B 75, 115335 

(2007) . 

A. Silva, Y. Oreg, and Y. Gefen, Phys. Rev. B 66, 195316 
(2002). 



5 



Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 W.-R. Lee, J. U. Kim, and H.-S. Sim, to be published. 

(1992). 



